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Abstract 

We discuss the dynamics of a spin coupled to a damped harmonic oscillator. This system can be mapped to a spin- 
boson model with a structured bath, i.e. the spectral function of the bath has a resonance peak. We diagonalize the 
model by means of infinitesimal unitary transformations {flow equations) , thereby decoupling the small quantum 
system from its environment, and calculate spin-spin correlation functions. 
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1. Introduction - Model 

Recently a new strategy for performing measure- 
ments on solid state (Josephson) qubits was proposed 
which uses the entanglement of the qubit with states 
of a damped oscillator [1], with this oscillator repre- 
senting the plasma resonance of the Josephson junc- 
tion. This system of a spin coupled to a damped har- 
monic oscillator (see Fig. 1) can be mapped to a stan- 
dard model for dissipative quantum systems, namely 
the spin-boson model [2]. Here the spectral function 
governing the dynamics of the spin has a resonance 
peak. Such structured baths were also discussed in con- 
nection with electron transfer processes [2] . We use the 
flow equation method introduced by Wegner [3] to an- 
alyze the system shown in Fig. 1, consisting of a two- 
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Fig. 1. A two-level system is coupled to a damped harmonic 
oscillator with frequency O. 



level system coupled to a harmonic oscillator fl, which 
is coupled to a bath of harmonic oscillators: 



H = — + nB^'B + giB'' + B)a, + ^ Cj^^blb 

k 

+ (fit + B) y K^ibi + + (B^ +Bfy^ 



with the spectral function J{ui) = k\5{u) — lik) = 
Tlo. This system can be mapped to a spin-boson 
model [2] 

n = -^'^x + \(J. ^ Mbl 6k) + ^ <-'k6l6k, (1) 



where the dynamics of the spin depends only on the 
spectral function J(cj) = — t^^k) given by 



J{u) 



with a - 



(2) 



2. Method - Results 

Using the flow equation technique we approximately 
diagonalize the Hamiltonian Ti. [Eq.(l)] by means of 
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infinitesimal unitary transformations. The continuous 
sequence of unitary transformations U (l) is labelled by 
a flow parameter I. Applying such a transformation 
to a given Hamiltonian, this Hamiltonian becomes a 
function of I: n{l) = U{l)nU\l). Here H{1 = 0) = H 
is the initial Hamiltonian and Ti.{l = oo) is the final 
diagonal Hamiltonian. Usually it is more convenient to 
work with a differential formulation 

^ = [r,il),nil)] with rj{l) = ^U-\l). (3) 

Using the flow equation approach one can decouple 
system and bath by diagonalizing Ti[l = 0) [4]: 



H{1 = oo) = -^^^ffx + y^^ujkblbk. 



(4) 



Here A^o is the renormalized tunneling frequency. For 
the generator of the flow we choose the Ansatz [4]: 

r, = J2 (i'^yAlft, + fct ) + a.c.k(6, - )) ^ (^^) 



(5) 



q,k 



with 7(tJk,a;q,Z) 



U)^ — UJq \u>h + A + A 



The flow equations for the effective Hamiltonian 
[Eq. (4)] then take the following form: 



^^^'^'^)=-2(a;-Afj(a;,0 



dl 



(6) 



+ 2AJ(a;,/) j dui'J{u!',l)I{uj,uj',l), 
= -A j dwJ{uj,l)-j-^. (7) 



The unitary flow diagonalizing the Hamiltonian gener- 
ates a flow for (7z(0 which takes the structure 



<T.(0 = h{r)<j^ + o-x ^Xk(0(&k + bi), 



(8) 



where h{l) and Xk(0 obey the differential equations 
Wk - A 



f = -A5^AkXk 



Wk + A' 



(9) 



^=AftAki^^ii-^+VxqAkAqA7(a;k,a;q,0. (10) 
al Wk + A ' 

q 

One can show that the function h{l) decays to zero as 
I ^ oo. Therefore the observable <Tz decays completely 
into bath operators [4] . 
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Fig. 2. (a) Diflcrcnl cflcctivc spectral functions ./(oj, / — 0) and 
(b) the corresponding C(cj) for OP — 0.06 and a — 0.15. The 
inset shows the term scheme of a two-level system coupled to 
a harmonic oscillator for the two limits Ao ^ and Ao ^ 

We integrated the flow equations numerically in or- 
der to calculate the Fourier transform, C(a;), of the 
spin-spin correlation function 



C{t) = i(o-z(i)<T.(0) + CTz(0)<T.(t)). 



(11) 



C{t) can be used to calculate dephasing and relaxation 

times for measurements on qubits [1]. Fig. 2(a) shows 
J(a', I = 0) and Fig. 2(b) C{lu) for different values of Q.. 
C{ijj) displays a double-peak structure, which can be 
understood from the term scheme shown in the inset. 
The arrows indicate the transitions responsible for the 
peaks in Ciuo). Additional structure of Ciuj) due to 
higher order transitions in the term scheme is not seen 
in Fig. 2. This is due to our Ansatz for ctz (0 [see Eq. (8)] , 
which does not include the corresponding higher order 
terms. However, we do not expect the additional peaks 
to have much weight, as the sum rule [4] for the total 
spectral weight is fulfilled with an error of less than 5% 
for all the plots in Fig. 2(b). We leave the extension of 
the Ansatz for crz(/) for future work. 
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